Relativistic-Newtonian correspondence of the zero-pressure but 
weakly nonlinear cosmology 

Hyerim Noh( a ) and Jai-chan Hwang( b ) 
( a ) Korea Astronomy and Space Science Institute, Daejon, Korea 
Department of Astronomy and Atmospheric Sciences, Kyungpook National University, Taegu, Korea 
E-mails: ^ hr@kasi.re.kr, ^ jchan@knu.ac.kr 
(February 7, 2008) 



in 
o 
o 

(N 



It is well known that couplings occur among the scalar-, vector-, and tensor-type perturbations of 
Friedmann world model in the second perturbational order. Here, we prove that, except for the 
gravitational wave contribution, the relativistic zero-pressure irrotational fluid perturbed to second 
order in a flat Friedmann background coincides exactly with the Newtonian result. Since we include 
the cosmological constant, our results are relevant to currently favoured cosmology. As we prove 
that the Newtonian hydrodynamic equations are valid in all cosmological scales to the second order, 
our result has an important practical implication that one can now use the large-scale Newtonian 
numerical simulation more reliably even as the simulation scale approaches and even goes beyond 
the horizon. That is, our discovery shows that, in the zero-pressure case, except for the gravitational 
wave contribution, there are no relativistic correction terms even near and beyond the horizon to 
the second-order perturbation. 
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I. INTRODUCTION 

Historically, the first proper cosmological analysis ap- 
peared only after the advent of Einstein's gravity theory 
in 1917 [1]: the only known preceding cosmologically rel- 
evant discussions can be found in Newton's correspon- 
dences to Bentley in 1692 [2]. The expanding world 
model and its linear structures were first studied in the 
context of Einstein's gravity in the classic studies by 
Friedmann in 1922 [3] and Lifshitz in 1946 [4], respec- 
tively. In an interesting sequence, the much simpler and, 
in hindsight, more intuitive Newtonian studies followed 
later by Milne in 1934 [5] and Bonnor in 1957 [6], re- 
spectively. According to Ellis "It is curious that it took 
so long for these dynamic models to be discovered af- 
ter the (more complex) general relativity models were 
known" [7]. This is particularly so, because in the case 
without pressure the Newtonian results coincide exactly 
with the previously derived relativistic ones for both the 
background world model and its first (linear) order per- 
turbations. It would be fair to point out, however, that 
the ordinarily known Newtonian cosmology (both for the 
Friedmann background and its linear perturbations) is 
not purely based on Newton's gravity, but is a guided one 
by Einstein's theory [8]. The zero-pressure system with 
the cosmological constant describes the current stage of 
our universe and its large scale structures in the linear 
stage remarkably well. 

Here we show that such a relativistic-Newtonian corre- 
spondence continues even to the weakly nonlinear order. 
As the observed large-scale structures show weakly non- 
linear processes our relativistic result has theoretical as 
well as practical significances to interprete and analyse 
such structures properly in relativistic level. As a conse- 
quence, our result implies that even to the weakly nonlin- 



ear (second perturbational) order the well known New- 
tonian equations can be used in all cosmological scales 
including the super-horizon scale. 
The known equations are 
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for the linear order perturbations [4,6]. The variable a(t) 
is the scale factor, and 5 = Sg/g with g and Sg the back- 
ground and perturbed parts of the density field, respec- 
tively; A is the cosmological constant. The "const." part 
is interpreted as the spatial curvature in Einstein's grav- 
ity [3,9] and the total energy in the Newton's gravity [5]. 
Although eq. (2) is also valid with general spatial curva- 
ture the relativistic-Newtonian correspondence is some- 
what ambiguous in the case with curvature [10]. There- 
fore, in the following we consider the flat background 
only. Equation (2) is valid even in the presence of the 
cosmological constant A. We will include A term in our 
analyses of the weakly nonlinear stage. The above two 
equations with vanishing spatial curvature describe re- 
markably well the current expanding stage of our uni- 
verse and its large-scale structures which are believed to 
be in the linear stage. In the small scale, however, the 
structures are apparently in nonlinear stage, and even in 
the large scale weakly nonlinear study is needed. Up un- 
til now, such a weakly nonlinear stage has been studied 
based on Newton's gravity only. 

The case with non-vanishing pressure cannot be han- 
dled in the Newtonian context, especially for the pertur- 
bation. In this work, we will show an additional contin- 
uation of relativistic-Newtonian correspondences in the 
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zero-pressure medium by showing that, except for the 
gravitational wave contribution, the relativistic second- 
order perturbation is described by the same set of equa- 
tions known in the Newtonian system. That is, except 
for the coupling with the gravitational waves, the Newto- 
nian equations coincide exactly with the relativistic ones 
even to the second order in perturbations. 

For relativistic perturbations, due to the covariance of 
field equations we have the freedom to fix the spacetime 
coordinates by choosing some of the metric or energy- 
momentum variables at our disposal: this is often called 
the gauge choice. The relativistic-Newtonian correspon- 
dence to the linear order was made by properly arrang- 
ing the equations using various gauge-invariant variables 
[11,10]. In the relativistic case 5 in eq. (2) is in fact 
5 in the temporal comoving gauge condition which also 
implies the temporal synchronous gauge in our zero- 
pressure case [4,12,11]. In this work we extend such cor- 
respondences to the second order. 

II. FULLY NONLINEAR EQUATIONS 

We may start from the completely nonlinear and co- 
variant equations; for a complete set of the covariant 
(1 + 3) equations see [14]. We consider a zero-pressure 
fluid with vanishing isotropic pressure and anisotropic 
stress, thus p = = n ab . In the energy frame we take 
q a = 0. Tildes indicate the covariant quantities, and the 
Greek and Latin indices indicate the space and spacetime 
indices, respectively. The momentum conservation gives 
vanishing acceleration vector a a to all orders; see eq. (27) 
of [13]. The energy conservation and the Raychaudhury 
equation (G^ — Gg part of Einstein's equation) give 

11 + = 0, (3) 

6 + + d ah ~o ah - u ab u ab + AirGf, - A = 0, (4) 

see eqs. (26,28) of [13]; 9 = u a . a is an expansion scalar 
with u a a fluid four-vector, a ab and u a b are the shear and 
the rotation tensors, respectively; if u a = we have no 
rotation of the fluid four- vector u a . We set c = 1. An 
overdot with tilde is a covariant derivative along the u a 
vector, e.g., fi = fL, a u a . By combining these equations 
we have 




- a ab a ab + u ab u ab - AnGfi + A = 0. 



(5) 

Equations (3-5) are valid to all orders, i.e., these equa- 
tions are fully nonlinear and still covariant. 

We consider the scalar- and tensor-type perturbations 
in the flat Friedmann background without pressure. We 



ignore the vector-type perturbation; the vector-type per- 
turbation (rotation) is supposed to be not important be- 
cause it always decays in the expanding phase even to 
the second order, see §VII.E of [13]. We will take the 
comoving gauge, and by ignoring the vector-type pertur- 
bations we have no rotation. In this case the four-vector 
becomes u a = 0, thus coincides with the normal four- 
vector h a . We lose no generality by imposing the gauge 
condition. In our case the energy-momentum tensor be- 
comes = — fi, and T" = = Tp where jl is the energy 
density We emphasize that as our comoving gauge con- 
dition fixes the gauge degree of freedom completely, all 
variables in this gauge condition are equivalcntly gauge 
invariant to the second order: this is in the sense that 
each of the variables has a unique corresponding gauge- 
invariant combination, see [13]. 

In [13] the equations are presented in the normal frame 
h a with n a = 0. The fluid four-vector u a in general 
differs from the normal four- vector h a . Only in the co- 
moving gauge without rotation the two frames coincide. 
Since the fluid quantities are defined in the fluid (u a ) 
frame, the zero-pressure condition should be imposed in 
u a frame. Thus, for the fluid quantities defined in the 
normal frame the physical zero-pressure condition implies 
vanishing pressures (both isotropic and anisotropic) only 
in the comoving gauge without rotation. In this nor- 
mal frame, the gauge transformation to the second order 
causes pressure terms to appear in other gauges, see [13]. 
In the energy frame, which takes vanishing flux q a = 
as the frame condition, the comoving gauge condition 
takes u a = for the fluid four-vector; here, we ignore 
the vector-type perturbation. Since u a — it coincides 
with the normal frame vector. Now, in the normal frame, 
which takes n a = as the frame condition, the comoving 
gauge condition without rotation implies q a = 0. Thus, 
as long as we take the comoving gauge without rotation, 
in either frame we have q~ a = and u a = = n a ; i.e., the 
fluid four-vector coincides with the normal four-vector. 



III. CORRESPONDENCE TO THE SECOND 
ORDER: A PROOF 

Now, we consider equations perturbed to the second 
order in the metric and matter variables. We introduce 

jl = IJ, + Sfi, 9 = 3- + 69, (6) 
a 

where \i and Sfi are the background and perturbed en- 
ergy density respectively, and 69 is the perturbed part 
of expansion scalar; we set 6 = 6/j,//i. We identify /j, = g 
to the background order, and 

5fi = 5g, 69 = ^-V ■ u, (7) 

to both the linear- and second-order perturbations. Now, 
to the second order, after some algebra using perturbed 
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order quantities presented in [13], the perturbed parts of 
eqs. (3,4) give (see the Appendix) 
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where cf^ is the transverse and tracefree tensor-type 
perturbation (the gravitational waves) introduced in eq. 
(15); the indices of are raised and lowered by 5 a p\ 
its contribution in cq. (9) comes from the shear term in 
eq. (4); see the Appendix. By combining these equations 
we have 
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which also follows from eq. (5). Equations (8-10) are our 
extension of eq. (2) to the second-order perturbations in 
Einstein's theory. We will show that, except for the grav- 
itational wave contribution, exactly the same equations 
also follow from Newton's theory. 

The presence of linear-order gravitational waves can 
generate the second-order scalar-type perturbation by 
generating the shear terms. The coupling between the 
scalar-type perturbation and the gravitational waves to 
the second order was noticed in the original study of the 
second-order perturbations by Tomita in 1967 [15]. Here, 
we note the behaviour of the gravitational waves in the 
linear regime. To the linear order the gravitational waves 
are described by the well known wave equation [4] 
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In the super-horizon scale the non-transient mode of C ( 

remains constant, thus C^l = 0, and in the sub-horizon 

scale, the oscillatory cfp redshifts away, thus Cfp oc a -1 . 

Notice that only time derivatives of C^l generate the 
scalar-type perturbation. Thus, we anticipate that the 
contribution of gravitational waves to the scalar- type per- 
turbation is not significant to the second order. The 
quadratic combinations of linear-order scalar-type per- 
turbation can also work as sources for the gravitational 
waves to the second-order [16]. 

In the Newtonian context, the mass conservation, the 
momentum conservation, and the Poisson's equation give 
[17] 
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V 2 <5$ = A-kGq5, 



(13) 
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where 5$ is the perturbed gravitational potential. We 
note that these equations are valid to fully nonlinear or- 
der. Equation (8) follows from eq. (12), and eq. (9) ig- 
noring the gravitational waves follows from eqs. (13,14). 
Thus, eq. (10) also naturally follows in Newton's the- 
ory [17]. In this Newtonian situation u is the perturbed 
velocity and 5 = Sg/g. This completes our proof of the 
relativistic-Newtonian correspondence to the second or- 
der. Although we successfully identified the relativistic 
density and velocity perturbation variables we do not 
have a relativistic variable which corresponds to the New- 
tonian potential <5<I> to the second order. This situation 
could be understood because the gravitational potential 
introduced in Poisson's equation reveals the action-at-a- 
distance nature and the static nature of Newton's gravity 
theory compared with Einstein's gravity. 

Notice that in the Newtonian context eqs. (8,9), thus 
eq. (10) as well, without the gravitational waves, are 
valid to fully nonlinear order. This has an important 
implication that any nonvanishing third and higher or- 
der correction terms in relativistic context should be re- 
garded as purely relativistic effects [18]. To our knowl- 
edge the Newtonian equations in eqs. (12-14) were first 
presented by Peebles in [17]. 



IV. DISCUSSIONS 

We have shown that to the second order, except for the 
gravitational wave contribution, the zero-pressure rela- 
tivistic cosmological perturbation equations can be ex- 
actly identified with the known equations in Newton's 
theory. As a consequence, to the second order, we iden- 
tified the correct relativistic variables which can be in- 
terpreted as density S/i and velocity 59 perturbations in 
eq. (7), and we showed that to that order the New- 
tonian hydrodynamic equations remain valid in all cos- 
mological scales including the super-horizon scale. Our 
results, showing the equivalence to the second order 
of the zero-pressure relativistic scalar-type perturbation 
and the Newtonian ones, may not be entirely surprising 
considering the Birkhoff's theorem [19]; for cosmology re- 
lated discussions of the theorem, see [17]. However, our 
results should not be so obvious either, because our sys- 
tem lacks any spatial symmetry contrary to the Birkhoff's 
theorem which is concerned with the spherically symmet- 
ric system. It might be as well that our relativistic results 
give relativistic correction terms appearing to the second 
order which become important as we approach and go be- 
yond the horizon scale. Our results show that there are 
no such correction terms appearing to the second order, 
and the correspondence is exact to that order. A comple- 
mentary result, showing the relativistic-Newtonian cor- 
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respondence in the Newtonian limit of post-Newtonian 
approach (^-expansion with ~ ^ < 1, thus valid 
far inside the horizon), can be found in [20]. In fact, 
the Newtonian hydrodynamic equations naturally appear 
in the zeroth post-Newtonian order of Einstein's gravity 
[21]; for the cosmological extension, see [22]. 

In a classic study of the cosmic microwave background 
radiation anisotropy in 1967 Sachs and Wolfe have men- 
tioned that "the linear perturbations are so surpris- 
ingly simple that a perturbation analysis accurate to sec- 
ond order may be feasible using the methods of Hawk- 
ing (1966)" [23]. Our proof of the exact relativistic- 
Newtonian correspondence to the second order could be 
regarded as one of such accurate results anticipated in 
[23] . Indeed, in this work we used the method of Hawk- 
ing which is the covariant (1+3) equations [24]; for other 
proofs see [16]. 

As we consider a flat background the ordinary Fourier 
analysis can be used to study the mode-couplings as in 
the Newtonian case [25]. Our equations include the cos- 
mological constant, thus compatible with current obser- 
vations of the large-scale structure and the cosmic mi- 
crowave background radiation which favour nearly flat 
Friedmann world model with non- vanishing A [26]. Our 
result may also have the following important practical 
cosmological implication. As we have proved that the 
Newtonian hydrodynamic equations are valid in all cos- 
mological scales to the second order, our result has an 
important comological implication that one can use the 
large-scale Newtonian numerical simulation more reliably 
in the general relativistic context even as the simulation 
scale approaches near (and goes beyond) horizon scale. 
The fluctuations near the horizon scale are supposed to 
be linear or weakly nonlinear; otherwise, it is difficult 
to imagine the presence of spatially homogeneous and 
isotropic background world model which is the basic as- 
sumption and the backbone of the modern cosmology. 

Since the Newtonian system is exact to the second or- 
der in nonlincarity, any non-vanishing third and higher 
order perturbation terms in the relativistic analysis can 
be regarded as the pure relativistic correction. Expand- 
ing the fully nonlinear equations in (3-5) to third and 
higher order will give the potential correction terms. For 
our recent work in the third-order perturbations, see 
[18]. In [18] we derive the nonvanishing third-order terms 
which are purely relativistic correction terms. We also 
show that these correction terms are independent of the 
horizon and are smaller than the second-order Newto- 
nian/relativistic terms by a factor 10~ 5 , thus negligible 
indeed. 



Acknowledgments 

We thank Dirk Puetzfeld for helpful comments. We 
wish to thank Lev Kofman and Ethan Vishniac for their 
comments on possible relation to the Birkhoff 's theorem. 



HN and JH were supported by grants No. R04-2003- 
10004-0 and No. R02-2003-000-10051-0, respectively, 
from the Basic Research Program of the Korea Science 
and Engineering Foundation. 



APPENDIX 

Since eqs. (8,9) allow us to conclude about the 
relativistic-Newtonian correspondence, in the following 
we will derive these equations in some detail using the 
basic quantities presented in [13]. As the metric we take 



ds 2 



-a 2 (1 + 2a) df] 2 - 2ax, a dr]dx° 



(l + 2p) 5 aP + 2C { *l 



dx a dx< 3 , 



(15) 



which follows from our convention in eqs. (49,175,178) 
of [13]. Here, a, \ an d <p are spacetime dependent 
perturbed-order variables; C^l is a transverse-tracefree 
tensor-type perturbation variable. In our metric we took 
the spatial C-gauge which removes the spatial gauge 
modes completely, thus all the remaining variables we are 
using arc spatially gauge- invariant, see §VI.B.2 of [13]; to 
the linear order, our notation coincides with Bardcen's 
in [27]. We work in the temporal comoving gauge which 
takes v = in [13]. The momentum conservation gives 
d a = which gives a = —^X' a X.a, see eq. (69) of [13]. 
In our comoving gauge condition, the four-vector in cq. 
(53) of [13], using eq. (175) in that paper, becomes 



u 



u° 



0; 



(l-2^-2C (t) ; 



(16) 



With these, we have 



// = fi.ou + jl^u a = [fj, (1 + 6)\ + -^/j,S^x a , 



1 3 H V u + 



(17) 



Our 59 is the same as — k in [13]. Using eqs. (55,57,70) 
of [13] we can show 



C dab 



x af> x,«p - 3 (Ax) 2 



+C {t)a "(^X, a0 + C^i). (18) 

Now, we have to relate \ to our notation. Apparently, 
we need \ onr y to the linear order. The G°-component 
of Einstein's equation in eq. (197) of [13] gives = 
— k = 59 = -V • u. As our u is of the potential type, 
i.e., of the form u = u, a , we have u = to the linear 
order. Using these, eqs. (3,4) give 
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( £ + 3-Vl + 6) + 5 + -V • u = --V • (Su) , (19) 
\fj, a J a a 

3- + AkGii - A + - V ■ | u + -u | + iirGuS 
a a \ a ) 

= "^V (u • Vu) - CW"' 3 g Uai/3 + . (20) 

To the background order, we have fi + 3(a/a)n — and 
3a/ a + A-kG^j, — A = 0; after an integration we recover eq. 
(1) with A; in this case the "const." is an integration con- 
stant which can be interpreted as the spatial curvature. 
The perturbed parts give eqs. (8,9). 
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